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Abstract 

We report first-principles investigations of the coloration of WO3 upon charge insertion, using 
sodium tungsten bronze (Na^WOs) as a model system. Our results explain well the systematic 
color change of Na x W03 from dark blue to violet, red-orange, and finally to golden-yellow as 
sodium concentration x increases from 0.3 to unity. Proper accounts for both the interband and 
the intraband contributions to the optical response are found to be very important for a detailed 
understanding of the coloration mechanism in this system. 
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Electro chromic materials exhibit reversible and persistent change of the optical 

properties, hence their color, upon applying an electrical pulse that injects both electrons 
and compensating ions into the material. Potential applications of electrochromic materials 

-jit— — a a 

Tungsten trioxide |3j, |4( is one of the most extensively studied electrochromic materials due 
to its superior coloration efficiency, short response time and reversibility Enhanced 
electrochromic properties in WO3 nanowires have been reported recently 5|, |6J. Despite 
much research effort, a first-principles theory for the coloration mechanism in this material 
has not emerged. Here we report density functional theory (DFT) 3] investigations of the 
coloration mechanism of WO3 upon charge insertion. Our results explains very well the 
systematic change in color S [9] of Na^WOa from dark blue to metallic golden-yellow 
with increasing sodium concentration. We find that proper accounts for the free-carriers 
contribution to the optical response are critical for a quantitative understanding of the 
coloration mechanism in this system. 

Undoped stoichiometric WO3 is a transparent semiconductor. Double injection of elec- 
trons and compensating ions (e.g., H + , Li + , or Na + ) by an electric pulse induces a dark blue 
color in the electrochromic WO3 film. Despite intense research over t he l ast three decades, 
the mechanism of electrochromism remains controversial 
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Interest- 



ingly, sodium tungsten bronze (Na^WOs) also shows various bluish colors at low sodium 
concentrations. Therefore, it is likely that both systems share the same fundamental col- 
oration mechanism 11], and the formation of tungsten bronze might be responsible for the 
electrochromism in WO3. Furthermore, the color of Na x W03 changes from dark blue to 
violet, red-orange, and finally to golden-yellow as sodium concentration x increases from 
0.3 to unity . A consistent theory should be able to explain the color of Na x W03 

for a wide range of Na concentration. Unfortunately, although the electronic and structural 
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properties of Na^WOs has been studied extensively 
quantitative understanding of its vivid color change with varying Na content is still lacking. 

Materials color from first principles — The apparent color of a material is intimately 
related to its optical absorption spectra. With the advance of modern electronic structure 
methods, various optical properties can be routinely calculated nowadays. However, to 
the best of our knowledge, there have been no reports of calculating the apparent color of 
materials from first-principles. The complicated physiology of human color perception makes 
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calculating the apparent color of a material rather difficult. Fortunately, the color space 23] 
defined by the International Commission of Illumination (CIE), the CIE 1931 XYZ color 
space, provides a mathematical foundation for a quantitative description of colors. 

The CIE XYZ color space defines three color matching functions, x(X), f/(A), and z(X) 
as shown in Fig. [TJ The XYZ tristimulus values for an object illuminated by a light source 
with a spectral power distribution 1(A) are then given by 



X = J I(X)R(X)x(X)dX; 
Y= JI(X)R(X)y(X)dX; 
Z= J I(X)R(X)z(X)dX, 



(1) 



where R(X) [or R{uj)} is the wavelength (frequency) dependent reflectivity of the object. The 
CIE XYZ color space is related to the well-known CIE RGB color space by the following 
linear transformation: 
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The reflective color and the relative brightness of a material can be constructed once we 
obtain the tristimulus values X, Y, and Z, or alternatively, R, G, and B. 

For an ideal semi-infinite crystal, the normal incidence reflectivity is related to the 
complex index of refraction n and the macroscopic complex dielectric function e(u) = 
Ei(uj) + i£2{oj) via 



R(u) 



'n(uj) + 1 



n(uj) 



£ 1/2 M 



(3) 



The imaginary part of the dielectric function E2{ui) due to interband transitions can be 
calculated using DFT-based first-principles electronic structure techniques: 
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where A is the polarization of light, v is the velocity operator, \ck) (\vk)) denotes unoccu- 
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pied conduction states (occupied valence states). The real part of the dielectric function 
is related to the imaginary part through the Kramers-Kronig relations. In undoped semi- 
conductors, only the interband transitions contribute to the optical response. In metals or 
doped semiconductors, however, the free-carriers contributions to the optical response are 
very important. These contributions can be conveniently modeled on the basis of the Drude 
theory: 



.intra. 



uj 1 = Aim/m* 



(5) 



cl>(7 2 + u 2 ) 

In the above equations, n, m*, and I/7 are, respectively, the density, optical effective mass, 
and relaxation time of the free carriers. Combining both the interband and intraband 
contributions thus gives the imaginary part of the total electronic dielectric function e* otal = 



,inter 



+ e 
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Electronic structure and optical properties of Na±,o WO3 — Tungsten trioxide comprises 
of corner- sharing WC>6 octahedra l|, |24|. Several metal ions can be incorporated into the 
structure to form tungsten bronzes. Among them sodium tungsten bronze (Na^WOs) is the 



most ex 
in color 



pensively studied system owing to its unusual optical properties and dramatic changes 
8|, |9] with increasing x. Na x W03 assumes a simple cubic structure for Na contents 
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x > 0.4 and undergoes a sequence of structural changes with decreasing x [25, 
phase can be synthesized under nonequilibrium conditions down to x ~ 0.2 
study covers Na concentration 0.3 < x < 1.0. For simplicity, we assume a cubic structure for 
all Na concentrations studied. Our calculations are based on DFT within the local density 



approximation (LDA) [28[. In order to avoid inconvenient (i.e^ 



arge) unit cells for systems 



with fractional x, we adopt the virtual crystal approximation [291. |30| in which fractional Na 
compositions are modeled by pseudoatoms with appropriate nuclear charges and valence 
electrons. We use the experimental lattice constants of Na^WOs, which can be conveniently 
expressed by the Vegard's law a(x) = 3.7845 + 0.0821a; 26]. The Brillouin zone summation 
in Eqn. [4] is carried out on a dense 32 x 32 x 32 uniform grid with an energy broadening of 
0.3 eV. 

Figure [2] (a) shows the band structure of Nai.oW03, the ending system of the series 
studied. As mentioned previously, Nai.oWOs displays a distinct golden-yellow color. The 
electronic band structure suggests, however, that the strong interband transitions mainly 
occur across an energy window of 4 ~ 7 eV. This falls well above the visible energy window 
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and cannot explain its golden- yellow appearance. There are also inter-subband transitions as 
indicated by circles in Fig. [2] (a). However, the energies of these transitions are very low (less 
than 1.0 eV) and cannot be responsible for the optical response in the visible range neither. 
Therefore, a proper account for the free carriers (intra-band) contributions is critical for a 
quantitative understanding of the optical properties of this system. 

Figure[2] (b) shows the calculated reflectivity spectra R(u) of Nai.oWOs, with and without 
including the free-carriers response. The unscreened plasma frequency in Eqn. [5] is calculated 
using the measured effective mass 15] . It is clear that the reflectivity in the visible range 
(1.5 ~ 3.0 eV) is dominated by the free carriers contribution (which is renormalized by the 
interband transitions as explained below). The spectra is fairly uniform and featureless for 
photon energy 1.5 < E < 3.0 eV if the free carriers plasma response is not included. As 
a result, Nai.oWOs would appear dark gray (more precisely, slightly bluish dark gray as 
shown in the inset of Fig. [2] (b) without taking into account the free-carriers contributions. 
These contributions significantly enhance the reflectivity below 2.5 eV and suppresses the 
reflectivity around 3.0 eV, resulting in the golden- yellow appearance of Nai.oWOs as shown in 
the inset of Fig.[2](b). The onset of the free-carriers enhanced reflectivity corresponds to the 
renormalized (screened) plasma frequency u p , defined by Re[e _1 ((D p )] = 0. The renormalized 
plasma frequency is related to the bare value (defined in Eqn. [5]) by u p « [u p / e™ ter (0)] 1 / 2 , 
which can be measured in electron energy loss spectroscopy experiments. Figure [2] (c) shows 
the calculated loss spectra of Nai.oWOs, which agree well with experiment [li|. Both the 
plasmon and interband excitations can be clearly identified. 

Optical properties and the color of Na x WO3 — The results for Nai.oWOs encourages us 
to carry out a systematic study of the optical properties of Na^WOs, with the hope that 
a quantitative understanding of dramatic color change of Na x WOs with varying x might 
be within the reach of today's first-principles electronic structure techniques. Figure [3] (c) 
shows the calculated reflectivity for Na^WOs (0.3 < x < 1.0). The onset of the free- 
carriers enhanced reflectivity shifts to shorter wavelengths as x increases. This is because 
the screened plasma frequency u p increases with increasing free-carrier concentration. These 



results agree very well with the measurements by Brown and Banks 
al. 



and by Goldner et 



311 ] except for an overall rescaling of the absolute magnitude of the reflectivity (likely 



due to sample quality and/or diffuse loss of light intensity). 

With these results in hand, we are now ready to calculate the XYZ and RGB tristimulus 
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values (denned in Eqn. [T] and [2]) of Na^WOs, and hence their colors. We assume that the 
illumination light source is uniform, i.e., 7(A) = constant in Eqn. [TJ Figure [3] (b) shows the 
calculated RGB components as functions of x. Note that the RGB values are normalized 
to unity for perfect reflectors [R{\) = 100%]. At low Na concentrations (0.3 < x < 0.4), 
the blue component dominates. This couples with the overall low reflectivity to give a 
dark bluish color for these systems as shown in Fig. [3] (c). In order to better illustrate the 
color, the RGB tristimulus values are rescaled to give an even brightness in Fig. [3] (c) . The 
calculated values for the relative brightness are also shown in the figure. The brightness 
values are normalized such that a perfect reflector will have a brightness of unity. 

As x increases beyond 0.6, both the red and green components are greatly enhanced, 
resulting in the red-orange to golden-yellow metallic appearance of these systems (shown in 
Fig. [3] (c). For intermediate Na concentration (0.5 < x < 0.6), red and blue components 
dominate; the systems thus appear violet to red-violet. Therefore, our theory is able to 
explain the systematic color change of Na^WOs from dark blue to metallic golden-yellow as 
x increases from 0.3 to 1.0. We mention that the description of color is somewhat subjective 
and it is not surprising that slightly different descriptions (e.g., purple v.s. violet) l], [si, [9J 
of the apparent color of Na^WOs are seen in literature. 

In summary, we have carried out first-principles calculations of the optical properties, 
and hence the color, of Na^WOs. The apparent color is constructed using the CIE 1931 
color matching functions and the calculated reflectivity spectra. Our results clearly ex- 
plains the systematic color change of Na x W03 as x varies and provide direct evidence that 
the coloration of WO3 upon charge insertion is a result of the subtle interplay between 
the interband and intraband optical response. The reflectivity is greatly enhanced below 
the screened plasma frequency, which explains the vivid color change of Na x W03 as Na 
concentration increases. In addition, we demonstrate that first-principles electronic struc- 
ture theory can predict the apparent color of crystalline materials and may help to design 
materials with desired colors for various applications. 
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FIG. 1: The CIE XYZ color matching functions. The CIE defined three color-matching functions, 
known as x(A), y(A), and z(X). These functions can be considered as the spectral sensitivity 
functions of three linear light detectors that give the CIE XYZ tristimulus values. 
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FIG. 2: Band structure (top), reflectivity (middle) and loss function of Nai.oWOs (bottom). Some 
strong interband transitions are indicated by shaded arrows. Strong inter-subband transitions are 
indicated by circles. The reflectivity spectra are calculated with and without the free-carriers 
contributions. Inset shows the colors of Nai.oW03 before and after taking into account the free- 
carriers response. 
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FIG. 3: Calculated reflectivity spectra (top), the RGB tristimulus values (middle), and the color 
of Na x W03 under uniform illumination (bottom) 
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